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1 Introduction 



This article presents the measurements of the differential inclusive J/ip production cross- 
section as a function of the J/ip transverse momentum, and of the fraction of J/ip mesons 
coming from the decay of a 6-hadron in pp collisions at a centre-of-mass energy of 2.76 TeV. 
The study is based on a sample corresponding to an integrated luminosity of 71nb _1 
collected in March 2011 with an average of one visible pp interaction per recorded event. 
The main goal of this short run was to provide a reference for the study of Pb-Pb 
interactions carried out at the same centre-of-mass energy per nucleon-nucleon collision. 

Studies of J/ip production have been performed by the LHC experiments using data 
taken at y/s = 7 TeV [l]44] as well as at lower energies [Hj. The data at y/s = 2.76 TeV 
provide an extra measurement to test theoretical models of J/ip production in hadron 
collisions and are also used to obtain a measurement of J/ip production from 6-hadron 
decays. 

The LHCb detector [6] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
(VELO) surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream. Charged hadrons 
are identified using two ring-imaging Cherenkov detectors. Photon, electron and hadron 
candidates are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a system which consists of five stations of alternating layers of iron and 
multiwire proportional chambers, with the exception of the centre of the first station, 
which uses triple-GEM detectors. 

For the data used in this analysis, the VELO, which consists of two retractable halves 
surrounding the interaction region, was positioned during collisions with its sensitive area 
at a minimum distance of 13 mm from the beam instead of the nominal 8 mm. This was 
necessary to provide a larger aperture for the beam at the lower centre-of-mass energy of 
2.76 TeV. 

The trigger (7) consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, which applies a full event reconstruction. 
Only the triggers used in this analysis are described here. At the hardware trigger level, 
a single muon candidate with larger than 0.8GeV/c is required. In the first stage of 
the software trigger a simplified event reconstruction is applied and one requires a [i + [i~ 
candidate with invariant mass greater than 2.7GeV/c 2 . In the second stage a full event 
reconstruction is performed and only events with a /i + yU~ pair with invariant mass within 
120MeV/c 2 of the known J/ip mass p] are retained. 



1 



2 Event selection 



The analysis strategy is based upon that described in Ref. (Tj. Candidate J/ip mesons 
are formed from pairs of opposite-sign charged particles reconstructed in the fiducial 
region 2 < rj < 5 by the full tracking system using algorithms adapted to the VELO at its 
displaced position. Each particle must have pr above 0.7 GeV/c and be identified as a muon. 
The two muons are required to originate from a common vertex, and only candidates with 
a y 2 probability of the vertex fit larger than 0.5% are kept. Events are selected in which 
at least one primary vertex is reconstructed from at least three VELO tracks, excluding 
the two signal muon tracks from the J/ip decay. A VELO track is required to have at least 
three hits on a straight line in the radial strips of the detector. 

The Monte Carlo samples used for this analysis are based on the Pythia 6.4 genera- 
is 



tor [9] configured with the parameters detailed in Ref. 10 . The EvtGen package 11 
used to generate hadron decays, in particular for J/ip and 6-hadrons. The interaction of the 
generated particles with the detector and its response are implemented using the Geant4 



toolkit [12] as described in Ref. 13 . Radiative corrections to the decay J/ip — > 



are generated using the Photos package 14 . The simulated position of the VELO 
corresponds to that in the data. 



3 Cross-section determination 

The differential cross-section for J/ip production in a px bin is given by 

d(T _ N (J/ip -» 

dp T £ x e tot x B (J/ip ->■ x Ap T ' 

where N (J/ip — >■ is the number of observed J/ip — » yU + /i~ signal decays in the given 

bin, e to t the J/ip detection efficiency per px bin (including both acceptance and trigger), 
C the integrated luminosity, B (J/ip — > = (5.93 ± 0.06) x 10~ 2 (8] the branching 

fraction of the J/ip — > [i + fi~ decay, and Apx the px bin size. 

The number of signal J/ip mesons per px bin is determined from an extended unbinned 
maximum likelihood fit to the invariant mass distribution of the reconstructed J/ip 
candidates in the interval 3.0 < M MAt < 3.2GeV/c 2 , where the signal is described by a 
Crystal Ball function |l5] and the combinatorial background by an exponential distribution. 
Figure [T] shows the J/ip invariant mass distribution together with the fit results for each 
Pt bin, where results for 7 < px < 12 GeV/c are merged in the last bin. 

There are two main sources that contribute to the inclusive J/ip sample. Those 
produced at the pp collision point, either directly or from the decay of a directly produced 
higher mass charmonium state, are called prompt J/ip . The second source, J/ip from b, 
are those produced in the decay of a 6-hadron. Their production is displaced from the pp 
collision point because of the relatively large b lifetime. The two sources are statistically 
separated using the measured J/ip pseudo-decaytime, defined as 

f = ( z Jrt ~ z pv) x M W / 2 x 

Pz 
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Figure 1: Dimuon mass distributions, with fit results superimposed, in bins of px- Results for 
Pt > 7GeV/c are merged in the last bin. 



where Zj/$ and zpy are the positions along the beam axis of the J/ip decay vertex and of 
the primary vertex refitted after removing the two muon tracks from the J/ip candidate; p z 
is the measured J/ip momentum in the beam direction and Mj/^ the known J/tp mass j8]. 
Given that 6-hadrons are not fully reconstructed, the J/tp momentum is used instead of the 
exact 6-hadron momentum and the t z variable provides a good estimate of the 6-hadron 
decaytime. 

The fraction of J/ip from b is determined from a simultaneous fit to the total pseudo- 
decaytime t z and invariant mass. Due to the small number of J/ip candidates, the 

fraction of J/ip from b is computed over the full px interval from to 12 GeV/c. The 
signal decaytime distribution is described by a delta function at t z = for the prompt 
J/ip component and an exponential decay function for the J/ip from b component. The 
function describing the t z distribution of the signal is therefore 

e T t 

/signa,l(t Z ; jfp, fb, 7b) = f p S(t z ) + 8{t z )f h , (3) 

Tb 

where 8(t z ) is the step function, f p and /b are the fractions of prompt J/ip and J/i/j from b 
in the sample, and 7~b the 6-hadron pseudo-lifetime. In the fit, r& is fixed to the value of 
1.52 ps, as obtained from simulation. The prompt and b components of the signal function 
are convolved with a triple-Gaussian resolution function 

a _{t z -n) 2 at _{t z -n) 2 1 _ a _ at (*z-m) 2 
f vcs (t z ; t i,a 1 ,a 2 ,a 3 ,P,P') = —^e ^T + -^e + P e ^ . (4) 

V27T(Ji V 27TCT 2 v27TO"3 

The parameter /i is the bias of the £ z measurement, and /3 and 0' the fractions of the 
first two Gaussian functions. The background consists of random combinations of muons 
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Figure 2: Distributions of the data with projections of the fit superimposed for (left) the dimuon 
invariant mass and (right) t z . The thick blue line is the total fit function and the hatched area 
the background component. 



from semi-leptonic b and c decays, which tend to produce positive t z values, as well as of 
mis-reconstructed tracks from decays in flight of kaons and pions, which contribute both 
to positive and negative t z values. The background t z distribution is parameterised with 
an empirical function based on the shape obtained from the J/ip mass sidebands. It is 
taken as the sum of a delta function and three exponential components, two for positive t z 
and one for negative t z . The exponential parameter, tl, is common to the larger positive 
and negative lifetime exponential components. The explicit form is 

_ tz \tz\ 

Wt z ) = (i - h - k) s(t z ) + e{t z )h — + f L e —^, (5) 

and is convolved with the same resolution function / res as the signal. 
The function used to describe the t z distribution is therefore 

/ -- 

f(tz] fp, fb, /i, O-i, 02, CT 3 , P, ft , T b ) = / p 8{t z ) + f h + (1 - f p - f h ) / bc kg(4y , , , 

V r h J (6) 

where all parameters except r b are freely varied. The total fit function is the sum of 
the products of the mass and t z fit functions for the signal and background. Figure [2] 
shows the distributions of the dimuon invariant mass and t z with the projections of the 
fit superimposed. The invariant mass resolution is 13.0 ± 0.3MeV/c 2 . The parameter /i 
describing the bias of the t z resolution function is 2.3 ± 2.0 fs and the RMS of the t z 
resolution function is 84 fs. As a measure of the fit quality, a x 2 is calculated using a 
binned event distribution. The resulting fit probability for the t z distribution is 90%. The 
fit gives a total yield of 3399 ± 65 J/ip signal decays. 
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p T (JA|/) [GeV/c] 

Figure 3: Acceptance (orange circles), reconstruction (upward- facing red triangles), trigger 
(downward- facing green triangles) and total (blue squares) J/tp efficiency, as a function of pt- 
The efficiencies are evaluated from a Monte Carlo simulation in which the J/tp is produced 
unpolarised. 

The fraction of signal J/ip coming from 6-hadron decays is measured to be F\> = j^j~ b — 
(6.7 ± 0.6)%. An absolute correction of 0.4% is applied based on simulation to take into 
account a bias produced by events in which 6-hadron decay products, other than the 
muons from the J/ip , are wrongly used to reconstruct the primary vertex. This leads to 
the result Fb = (7.1 ± 0.6)% where the uncertainty is only statistical. 

A simulated sample of inclusive, unpolarised J/ip mesons is used to estimate the 
geometrical acceptance in each pr bin. The reconstruction efficiency, which combines 
the J/ip meson detection, reconstruction and selection efficiencies, is also computed from 
simulation as a function of p^ and is corrected to account for the difference observed in the 
tracking efficiency between data and simulation at -^s = 7 TeV. This correction is about 
1%. The efficiency of the hardware trigger is determined directly from data using a large 
inclusive J/ip sample at a/s = 7 TeV triggered and selected with the same requirements as 
those used in this analysis: the efficiency is calculated in small bins of the J/ip transverse 
momentum and rapidity and weighted according to the p^ and y distributions as given 
by the simulation at 2.76 TeV. The efficiency of the software trigger, which makes use of 
the VELO information, is determined from simulation since the data at ^/s = 7 TeV were 
taken with the VELO in the closed position. The total efficiency, calculated as the product 
of acceptance, reconstruction and trigger efficiencies, and its components are displayed in 
Fig. [3] as a function of px- A non-zero polarisation of the J/0 at production can affect 
the total efficiency [I]. The results quoted in this article assume that the J/ip mesons are 
produced unpolarised. 
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4 Luminosity determination 



To determine the integrated luminosity, an effective interaction rate is continuously 
measured during data taking and an absolute calibration is performed with a dedicated van 



der Meer (VDM) scan 16 . The strategy is similar to that developed for the y/s = 7TeV 
running [l~7] . 

The VDM method exploits the ability to move the beams in both transverse coordinates 
with high precision and thus to scan the colliding beams with respect to each other. The 
limiting systematic uncertainty affecting the VDM measurement arises from the knowledge 
of the number of protons in the colliding bunch pairs. These are measured with two 



types of beam current transformers installed in the LHC 18 . The DCCT (DC Current 



Transformer) measures the total beam current, and is thus used to constrain the total 



number of particles. The uncertainty associated with the DCCT calibration is 2.7% 19 
The other transformer, the FBCT (Fast Beam Current Transformer) is used to measure 
the relative charges of the individual bunches. The uncertainty in its offset and linearity 
contributes a 0.9% uncertainty to the overall luminosity [19]. A small fraction of protons 
in the LHC may be captured outside the nominally filled bunch slots. This contribution, 
which needs to be subtracted from the DCCT measurement, is estimated to be 2.5% from 
the number of beam-gas events in nominally empty bunch crossings. Due to the small 
number of such events and uncertainties in the trigger efficiency, the subtraction introduces 
a cross-section uncertainty of 0.4%. The uncertainty in the length-scale calibration, which 
affects the beam separation values, contributes 1% to the systematic uncertainty in the 
luminosity. Finally, a 2.1% uncertainty is assigned to account for a non-reproducibility of 
the VDM results observed when performing similar luminosity calibration measurements 



Table 1: Relative systematic uncertainties on the luminosity (%). 



Uncertainty on relative normalisation 

Counter stability 0.5 

/i variation among bunch crossings 0.5 

Uncertainty on absolute normalisation 

Statistical error of the VDM scan 0.2 

Total beam current 2.7 

Individual bunch population 0.9 

Protons outside nominal bunches 0.4 

Length scale calibration 1.0 

Non-reproducibility in similar scans at \/s = 7TeV 2.1 

Total uncertainty 3.8 
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at yfs = 7TeV, as decribed in Ref. fl7) . 

The integrated luminosity for the runs considered in this analysis is measured to be 
70.6±2.7nb _1 . A summary of the contributions to the overall luminosity uncertainty is 
provided in Table [1} The uncertainties are uncorrelated and therefore added in quadrature. 

5 Systematic uncertainties 

The different contributions to the systematic uncertainty affecting the cross-section mea- 
surement are summarised in Table [2j Correction factors estimated directly from data to 
take into account residual differences between simulation and data are also detailed. 

The influence of the choice of the fit function used to describe the shape of the dimuon 
mass distribution is estimated by fitting the J/ip invariant mass distribution with the sum 
of two Crystal Ball functions. The relative difference of 2.2% in the number of signal 
events is taken as systematic uncertainty. 

A fraction of J/ip events have a lower mass because of the radiative tail. Based on 
Monte Carlo studies, 5% of the J/ip signal is estimated to be outside the analysis mass 
window (M w < 3.0 GeV/c 2 ) and not counted as signal. The fitted signal yields are therefore 
corrected, and an uncertainty of 1% is assigned to the cross-section measurement based on 
a comparison between the radiative tail observed in data and simulation. 

To cross-check and assign a systematic uncertainty to the Monte Carlo determination of 
the muon identification efficiency, the single track muon identification efficiency is measured 
on data using a tag-and-probe method. This method reconstructs J/ip candidates in which 
one muon is identified by the muon system ( "tag" ) and the other one ( "probe" ) is identified 
by selecting a track with a minimum-ionising energy deposition in the calorimeters. The 
absolute muon identification efficiency is then evaluated on the probe muon, as a function 
of the muon momentum and found to be larger than 95%. The ratio of the muon 
identification efficiency measured in data to that obtained in the simulation is convolved 
with the momentum distribution of muons from J/ip to obtain an efficiency correction. 
This factor is found to be f .024 ± 0.011 and is consistent with being constant over the 
full J/ip transverse momentum and rapidity range; the error on the correction factor is 
included systematic uncertainty. 

Studies at \fs = 7TeV have shown that the Monte Carlo simulation reproduces the 
determination from data of the efficiency to reconstruct the two muon tracks from the J/ip 
decay within 0.8% to 1.1%, depending on the J/ip transverse momentum. This difference 
is taken as a systematic uncertainty. An additional uncertainty of 1% per track is assigned 
to cover differences in the efficiency of the track x 2 /ndf cut between data and simulation. 
Similarly, for the selection based on the J/ip vertex y 2 probability, a difference below 0.3% 
is measured between the cut efficiency computed in data and simulation, which is assigned 
as systematic uncertainty. To take into account the model dependence of the simulation 
in the efficiency calculation, the main parameters of the Pythia 6.4 generator related to 
prompt J/ip production were varied. These parameters define the minimum cut-offs for 
regularising the cross-section. A 4.5% effect on the total efficiency was observed. 
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Table 2: Relative systematic uncertainties on the cross-section results and on the fraction of 
J/ip mesons from 6-hadron decay (%). 



Source 


Systematic uncertainty 


Correlated between bins 




Mass fits 


2.2 


Radiative tail 


1.0 


Muon identification 


1.1 


Tracking efficiency 


0.8 to 1.1 


Track % 2 


2.0 


Vertexing 


0.3 


Model dependence 


4.5 




1.0 


Luminosity 


3.8 


Uncorrelated between bins 


Trigger 


1.6 to 7.7 


Applied only to J/ip from b fraction 


t z fit 


10.0 


Applied only to cr(pp — > bbX) 


B(b^ J/ipX) 


8.6 



The hardware trigger efficiency is determined using a sample of events at y/s = 7 TeV 
that would still be triggered if the J/0 candidate were removed. The software trigger 
efficiency is obtained from the simulation. Its uncertainty is evaluated by comparing true 
and measured trigger efficiency using a trigger-unbiased sample of simulated J/ip events. 

Uncertainties related to the t z fit procedure are taken into account by varying the 
slope of the exponential function of the J/ip from b component by its uncertainty in 
the simulation (2%). The resulting 10% variation of the number of J/ip from b is used 
as a systematic uncertainty that affects the measurement of F^. The influence of the 
background parametrisation was studied by varying the number of exponential functions 
in Eq. (JsJ> and found to be negligible. Furthermore, an uncertainty of 8.6% on the average 
branching fraction of b decays to a final state containing a J/ip meson contributes to the 
uncertainty on the extrapolation to the total bb cross-section. 
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Table 3: Differential cross-section da /dpi at y/s = 2.76 TeV for inclusive J/ip production in bins 
of pr- The rapidity range covered is 2.0 < y < 4.5. The first uncertainty is statistical and the 
second is systematic. 



p T (GeV/c) 


dcr/dpr [nb/(GeV/c)] 


0-1 


1270 ±60 ±130 


1-2 


1780 ± 70 ± 160 


2-3 


1290 ± 50 ± 90 


3-4 


700±40± 50 


4-5 


313±22± 24 


5-6 


142±13± 10 


6-7 


61 ± 8± 4 


7-12 


14 ± 2± 1 



6 Results 

The measured differential cross-section for inclusive J/ip production as a function of p^, 
after all corrections and assuming no polarisation, is given in Table [3] and displayed in 
Fig. |4j The integrated inclusive cross-section for J/ip production in the defined fiducial 
region is 

a (J/ip , p T < 12 GeV/c, 2.0 < y < 4.5) = 5.6 ± 0.1 ± 0.4 \ib. 

The first uncertainty is statistical and the second systematic. Studies indicate that this 
result could change by up to 20% assuming fully longitudinal or fully transverse J/ip 
polarisation |lj. The fraction of J/ip from b is measured to be 

F b = (7.1 ±0.6 ±0.7)% 

in the same acceptance range, < 12GeV/c and 2.0 < y < 4.5. 
From the above results, one can deduce 

a (J/ip from b, p T < 12 GeV/c, 2.0 < y < 4.5) = 400 ± 35 ± 49 nb, 

in good agreement with the theoretical prediction of 370^^° n b? based on NLO calculations 



described in Ref. 20 . In addition, the total bb production cross-section is computed as 



a (J/iP, p T < 12 GeV/c, 2.0 < y < 4.5) x F h 
a( PP -+ bbX) = a 47r 2B (b^J/ipX) ' (?) 

where the factor = 6.3 is an extrapolation factor of the cross-section from 

the measured to the full kinematic region. This factor is obtained using the sim- 
ulation software described previously. The inclusive b — > J/ipX branching frac- 
tion is B(b — >■ J/tpX) = (1.16 ± 0.10)% l8J. The resulting total bb cross-section is 
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Figure 4: Differential production cross-section for inclusive J/ip production in the rapidity range 
2.0 < y < 4.5 as a function of pr- The vertical error bars are the quadratic sums of the statistical 
and systematic uncertainties. 



a(pp — > bbX) = 110 ± 9 ± 16 |J.b. No systematic uncertainty has been included for the 
extrapolation factor estimated from the simulation. The value of the extrapolation 



factor given by NLO calculations is 6.1 20 . 



7 Conclusions 

The differential cross-section for inclusive J/ip production is measured as a function of the 
J/ip transverse momentum in the forward region, 2.0 < y < 4.5. The analysis is based on 
a data sample corresponding to an integrated luminosity of 71 nb _1 collected by the LHCb 
experiment at the Large Hadron Collider at a centre-of-mass energy of \/s = 2.76 TeV. 
The results obtained are in good agreement with earlier measurements of the inclusive 
J/ip production cross-section in pp collisions at the same centre-of-mass energy, performed 



by ALICE in the region 2.5 < y < 4.0 21 . A first measurement of the production of J/ip 



from 6-hadron decays at 2.76 TeV is also obtained. 
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